Enhancing agro-biodiversity in chicken: a sensory comparison of broths from German local chicken breeds and their crossbreeds.

The poultry meat production landscape has undergone a reduction in chicken breeds, resulting in a reliance on a limited number of varieties. Motivated by the goal of promoting sustainable chicken production and enhancing agro-biodiversity, this study pioneers a comparison between local chicken breeds (LB) and their crossbreeds (CB) with modern hybrid lines. Serving as an initial exploration within a larger project, this research acts as a prelude to a comprehensive investigation, aiming to complement the human sensory assessment of product quality. Study I assessed chicken broths prepared from 3 German LBs Bielefelder (BIE), Altsteirer (ALT), and Ramelsloher (RAM) utilizing a factorial 3 × 2 × 2 design that incorporated variations in salt content (unsalted/salted) and cooking time (1 h/3 h). The sensory profiles of the LB broths were largely similar, except for BIE, which exhibited a higher skin odor intensity. Both, increased salt content and longer cooking time intensify sensory perception on most attributes. In study II, BIE was compared with 6 CBs, with variations in salt content and cooking time (6 × 2 × 2 + 1 × 2 factorial design). BIE demonstrated higher sensory intensities than the CBs. Those were comparable, with no clear advantages or disadvantages identified from a sensory standpoint. These findings support that crossbreeding with commercial lines is not associated with changes in the sensory profile. It thus represents a strategy for improving the economic viability of local chicken breeds in order to preserve their valuable agro-biodiversity. The provided protocol for evaluating chicken broth from LBs or their CBs aims to offer researchers a standardized foundation for sensory assessments in chicken broth studies.

Mendelian randomization analysis of 34,497 German Holstein cows to infer causal associations between milk production and health traits.

BACKGROUND: Claw diseases and mastitis represent the most important health issues in dairy cattle with a frequently mentioned connection to milk production. Although many studies have aimed at investigating this connection in more detail by estimating genetic correlations, they do not provide information about causality. An alternative is to carry out Mendelian randomization (MR) studies using genetic variants to investigate the effect of an exposure on an outcome trait mediated by genetic variants. No study has yet investigated the causal association of milk yield (MY) with health traits in dairy cattle. Hence, we performed a MR analysis of MY and seven health traits using imputed whole-genome sequence data from 34,497 German Holstein cows. We applied a method that uses summary statistics and removes horizontal pleiotropic variants (having an effect on both traits), which improves the power and unbiasedness of MR studies. In addition, genetic correlations between MY and each health trait were estimated to compare them with the estimates of causal effects that we expected. RESULTS: All genetic correlations between MY and each health trait were negative, ranging from - 0.303 (mastitis) to - 0.019 (digital dermatitis), which indicates a reduced health status as MY increases. The only non-significant correlation was between MY and digital dermatitis. In addition, each causal association was negative, ranging from - 0.131 (mastitis) to - 0.034 (laminitis), but the number of significant associations was reduced to five nominal and two experiment-wide significant results. The latter were between MY and mastitis and between MY and digital phlegmon. Horizontal pleiotropic variants were identified for mastitis, digital dermatitis and digital phlegmon. They were located within or nearby variants that were previously reported to have a horizontal pleiotropic effect, e.g., on milk production and somatic cell count. CONCLUSIONS: Our results confirm the known negative genetic connection between health traits and MY in dairy cattle. In addition, they provide new information about causality, which for example points to the negative energy balance mediating the connection between these traits. This knowledge helps to better understand whether the negative genetic correlation is based on pleiotropy, linkage between causal variants for both trait complexes, or indeed on a causal association.

Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows.

BACKGROUND: Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM. RESULTS: In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect. CONCLUSION: Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.

Genomic evidence for the suitability of Göttingen Minipigs with a rare seizure phenotype as a model for human epilepsy.

Epilepsy is a complex genetic disorder that affects about 2% of the global population. Although the frequency and severity of epileptic seizures can be reduced by a range of pharmacological interventions, there are no disease-modifying treatments for epilepsy. The development of new and more effective drugs is hindered by a lack of suitable animal models. Available rodent models may not recapitulate all key aspects of the disease. Spontaneous epileptic convulsions were observed in few Göttingen Minipigs (GMPs), which may provide a valuable alternative animal model for the characterisation of epilepsy-type diseases and for testing new treatments. We have characterised affected GMPs at the genome level and have taken advantage of primary fibroblast cultures to validate the functional impact of fixed genetic variants on the transcriptome level. We found numerous genes connected to calcium metabolism that have not been associated with epilepsy before, such as ADORA2B, CAMK1D, ITPKB, MCOLN2, MYLK, NFATC3, PDGFD, and PHKB. Our results have identified two transcription factor genes, EGR3 and HOXB6, as potential key regulators of CACNA1H, which was previously linked to epilepsy-type disorders in humans. Our findings provide the first set of conclusive results to support the use of affected subsets of GMPs as an alternative and more reliable model system to study human epilepsy. Further neurological and pharmacological validation of the suitability of GMPs as an epilepsy model is therefore warranted.

Simulation of dual-purpose chicken breeding programs implementing gene editing.

BACKGROUND: In spite of being controversial and raising ethical concerns, the application of gene editing is more likely to be accepted when it contributes to improving animal welfare. One of the animal welfare and ethical issues in chicken breeding is chick culling, the killing of the male layer chicks after hatching due to the poor fattening performance. Although establishing dual-purpose chicken lines could solve this problem, unfavorable genetic correlations between egg and meat production traits hindered their competitiveness. Although it is also controversial in ethical terms, gene editing may accelerate genetic progress in dual-purpose chicken and alleviate the ethical concerns from chick culling. RESULTS: The simulation compared the utility improvement in dual-purpose use under two breeding schemes: one consisting in the improvement of the laying hens, and the second in the improvement of a synthetic line obtained from a layer broiler cross. In each breeding scheme, the breeding programs were simulated with and without gene editing. Polygenic breeding values and 500 simulated quantitative trait loci (QTL) with different levels of pleiotropy caused negative correlations between egg production, meat production, and overall health. The results of the simulation demonstrated that genetic gain could be accelerated by at most 81% for several generations if gene editing was used. The actual increase in genetic gain depended on the number of single nucleotide polymorphisms (SNPs) being edited per animal. The rate of genetic improvement became equal in scenarios with and without gene editing after 20 generations. This is because the remaining segregating QTL had small effects and their edition would have negative overall health effects from potential off-target edits. Although gene editing can improve genetic gain in quantitative traits, it can only be recommended as long as QTL with reasonable effect sizes are segregating and detectable. CONCLUSIONS: This simulation demonstrates the potential of gene editing to accelerate the simultaneous improvement of negatively correlated traits. When the risk of negative consequences from gene editing persists, the number of SNPs to be edited should be chosen carefully to obtain the optimal genetic gain.

Defining valid breeding goals for animal breeds.

BACKGROUND: The objective of any valid breeding program is to increase the suitability of a breed for its future purposes. The approach most often followed in animal breeding for optimizing breeding goals assumes that the sole desire of the owners is profit maximization. As this assumption is often violated, a generalized approach is needed that does not rely on this assumption. RESULTS: The generalized approach is based on the niche concept. The niche of a breed is a set of environments in which a small population of the breed would have a positive population growth rate. Its growth rate depends on demand from prospective consumers and supply from producers. The approach involves defining the niche that is envisaged for the breed and identifying the trait optima that maximize the breed's adaptation to its envisaged niche within the set of permissible breeding goals. The set of permissible breeding goals is the set of all potential breeding goals that are compatible with animal welfare and could be reached within the planning horizon of the breeding program. In general, the breed's adaptation depends on the satisfaction of the producers with the animals and on the satisfaction of the consumers with the products produced by the animals. When consumers buy live animals, then the breed needs to adapt to both the environments provided by the producers, and the environments provided by the consumers. The profit function is replaced by a more general adaptedness function that measures the breed's adaptation to its envisaged niche. CONCLUSIONS: The proposed approach coincides with the traditional approach if the producers have the sole desire to maximize their income, and if consumer preferences are well reflected by the product prices. If these assumptions are not met, then the traditional approach to breeding goal optimization is unlikely to result in a valid breeding goal. Using the example of companion breeds, this paper shows that the proposed approach has the potential to fill the gap.

In vitro modelling of the influence of alternative feeds (Hermetia illucens, Arthrospira platensis) on the resistance of different rainbow trout populations (Oncorhynchus mykiss) against the viral haemorrhagic septicaemia virus and Yersinia ruckeri.

Replacing fishmeal, a finite resource with high market demand, in the diet of carnivorous rainbow trout with proteins from alternative sources may be a challenge for these fish. Therefore, this study investigated whether replacing fishmeal with protein derived from Hermetia illucens or Arthrospira platensis could promote disease susceptibility in local trout populations with different growth performance. This was assessed in vitro by measuring susceptibility to infection with the viral haemorrhagic septicaemia virus (VHSV) or the bacterium Yersinia ruckeri. Analysis of fin tissue explants and primary cell cultures from scales from the three trout populations infected in vitro with VHSV and gill explants infected with Y. ruckeri showed no significant differences in virus replication or bacterial counts. Evaluation of the virucidal or bactericidal effect of skin mucus showed a significant reduction in viral load and bacterial count for all samples with mucus addition, but no significant difference was observed between the experimental groups. This study documents no apparent impairment of innate immune mechanisms in the skin and gills of trout after feeding a diet replacing fishmeal with Arthrospira or Hermetia proteins. This underlines the potential of these alternative protein sources for the further development of sustainable trout aquaculture.

Blood transcriptome analysis in a buck-ewe hybrid points towards an nuclear factor-kappa B lymphoproliferative autoimmune disorder.

Mammal hybridization is a speciation mechanism and an evolutionary driver. Goat-sheep, especially buck-ewe hybrids, are very rare with only one case reported in 2016, which is the subject of the work presented here. Blood transcriptome analysis revealed that the hybrid largely deviated from imprinting schemes previously described in sheep and other mammals. Furthermore, transcriptome regulation seems to differ from the parent transcriptomes, which is most likely a product of partially incompatible imprinting mechanisms from two closely related species. To gain a deeper understanding of hybridization in mammals we re-analyzed the RNA sequencing data of the buck-ewe hybrid and its parents. We found parent-of-origin-specific expression of genes that functionally clustered, which we explain with the Dobzhansky-Muller incompatibility (DMI) model. According to the DMI model, proteins which interact have a high probability of being barrier loci and hence are prone to monoallelic expression. We discovered enrichment of genes uniquely expressed by the buck-ewe hybrid, which implicate that it suffered from an NF-κB lymphoproliferative autoimmune disorder. Similar findings were reported in the F1 generation of hybrid mice. We propose that hybridization of two related species may lead to an autoimmune phenotype, due to immunoglobulin incompatibilities and incomplete silencing of barrier loci.

Single-step SNPBLUP evaluation in six German beef cattle breeds.

The implementation of genomic selection for six German beef cattle populations was evaluated. Although the multiple-step implementation of genomic selection is the status quo in most national dairy cattle evaluations, the breeding structure of German beef cattle, coupled with the shortcoming and complexity of the multiple-step method, makes single step a more attractive option to implement genomic selection in German beef cattle populations. Our objective was to develop a national beef cattle single-step genomic evaluation in five economically important traits in six German beef cattle populations and investigate its impact on the accuracy and bias of genomic evaluations relative to the current pedigree-based evaluation. Across the six breeds in our study, 461,929 phenotyped and 14,321 genotyped animals were evaluated with a multi-trait single-step model. To validate the single-step model, phenotype data in the last 2 years were removed in a forward validation study. For the conventional and single-step approaches, the genomic estimated breeding values of validation animals and other animals were compared between the truncated and the full evaluations. The correlation of the GEBVs between the full and truncated evaluations in the validation animals was slightly higher in the single-step evaluation. The regression of the full GEBVs on truncated GEBVs was close to the optimal value of 1 for both the pedigree-based and the single-step evaluations. The SNP effect estimates from the truncated evaluation were highly correlated with those from the full evaluation, with values ranging from 0.79 to 0.94. The correlation of the SNP effect was influenced by the number of genotyped animals shared between the full and truncated evaluations. The regression coefficients of the SNP effect of the full evaluation on the truncated evaluation were all close to the expected value of 1, indicating unbiased estimates of the SNP markers for the production traits. The Manhattan plot of the SNP effect estimates identified chromosomal regions harbouring major genes for muscling and body weight in breeds of French origin. Based on the regression intercept and slope of the GEBVs of validation animals, the single-step evaluation was neither inflated nor deflated across the six breeds. Overall, the single-step model resulted in a more accurate and stable evaluation. However, due to the small number of genotyped individuals, the single-step method only provided slightly better results when compared to the pedigree-based method.

Disturbed circadian rhythm of locomotor activity of pullets is related to feather pecking in laying hens.

Various aspects of activity, such as spontaneous activity, explorative activity, activity in open-field tests, and hyperactivity syndrome have been explored as causal factors of feather pecking in laying hens, with no clear results. In all previous studies, mean values of activity over different time intervals were used as criteria. Incidental observation of alternated oviposition time in lines selected for high (HFP) and low feather pecking (LFP), supported by a recent study which showed differentially expressed genes related to the circadian clock in the same lines, led to the hypothesis that feather pecking may be related to a disturbed diurnal activity rhythm. Hence activity recordings of a previous generation of these lines have been reanalyzed. Data sets of a total of 682 pullets of 3 subsequent hatches of HFP, LFP, and an unselected control line (CONTR) were used. Locomotor activity was recorded in pullets housed in groups of mixed lines in a deep litter pen on 7 consecutive 13-h light phases, using a radio-frequency identification antenna system. The number of approaches to the antenna system was recorded as a measure of locomotor activity and analyzed using a generalized linear mixed model including hatch, line, time of day and the interactions of hatch × time of day and line × time of day as fixed effects. Significant effects were found for time and the interaction line × time of day but not for line. All lines showed a bimodal pattern of diurnal activity. The peak activity of the HFP in the morning was lower than that of the LFP and CONTR. In the afternoon peak all lines differed with the highest mean in the LFP followed by CONTR and HFP. The present results provide support for the hypothesis that a disturbed circadian clock plays a role in the development of feather pecking.

Genomic dominance variance analysis of health and milk production traits in German Holstein cattle.

Genomic analyses commonly explore the additive genetic variance of traits. The non-additive variance, however, is usually small but often significant in dairy cattle. This study aimed at dissecting the genetic variance of eight health traits that recently entered the total merit index in Germany and the somatic cell score (SCS), as well as four milk production traits by analysing additive and dominance variance components. The heritabilities were low for all health traits (between 0.033 for mastitis and 0.099 for SCS), and moderate for the milk production traits (between 0.261 for milk energy yield and 0.351 for milk yield). For all traits, the contribution of dominance variance to the phenotypic variance was low, varying between 0.018 for ovarian cysts and 0.078 for milk yield. Inbreeding depression, inferred from the SNP-based observed homozygosity, was significant only for the milk production traits. The contribution of dominance variance to the genetic variance was larger for the health traits, ranging from 0.233 for ovarian cysts to 0.551 for mastitis, encouraging further studies that aim at discovering QTLs based on their additive and dominance effects.

A genomic assessment of the correlation between milk production traits and claw and udder health traits in Holstein dairy cattle.

Claw diseases and mastitis represent the most important disease traits in dairy cattle with increasing incidences and a frequently mentioned connection to milk yield. Yet, many studies aimed to detect the genetic background of both trait complexes via fine-mapping of quantitative trait loci. However, little is known about genomic regions that simultaneously affect milk production and disease traits. For this purpose, several tools to detect local genetic correlations have been developed. In this study, we attempted a detailed analysis of milk production and disease traits as well as their interrelationship using a sample of 34,497 50K genotyped German Holstein cows with milk production and claw and udder disease traits records. We performed a pedigree-based quantitative genetic analysis to estimate heritabilities and genetic correlations. Additionally, we generated GWAS summary statistics, paying special attention to genomic inflation, and used these data to identify shared genomic regions, which affect various trait combinations. The heritability on the liability scale of the disease traits was low, between 0.02 for laminitis and 0.19 for interdigital hyperplasia. The heritabilities for milk production traits were higher (between 0.27 for milk energy yield and 0.48 for fat-protein ratio). Global genetic correlations indicate the shared genetic effect between milk production and disease traits on a whole genome level. Most of these estimates were not significantly different from zero, only mastitis showed a positive one to milk (0.18) and milk energy yield (0.13), as well as a negative one to fat-protein ratio (-0.07). The genomic analysis revealed significant SNPs for milk production traits that were enriched on Bos taurus autosome 5, 6, and 14. For digital dermatitis, we found significant hits, predominantly on Bos taurus autosome 5, 10, 22, and 23, whereas we did not find significantly trait-associated SNPs for the other disease traits. Our results confirm the known genetic background of disease and milk production traits. We further detected 13 regions that harbor strong concordant effects on a trait combination of milk production and disease traits. This detailed investigation of genetic correlations reveals additional knowledge about the localization of regions with shared genetic effects on these trait complexes, which in turn enables a better understanding of the underlying biological pathways and putatively the utilization for a more precise design of breeding schemes.

Structural variation and eQTL analysis in two experimental populations of chickens divergently selected for feather-pecking behavior.

Feather pecking (FP) is a damaging nonaggressive behavior in laying hens with a heritable component. Its occurrence has been linked to the immune system, the circadian clock, and foraging behavior. Furthermore, dysregulation of miRNA biogenesis, disturbance of the gamma-aminobutyric acid (GABAergic) system, as well as neurodevelopmental deficiencies are currently under debate as factors influencing the propensity for FP behavior. Past studies, which focused on the dissection of the genetic factors involved in FP, relied on single nucleotide polymorphisms (SNPs) and short insertions and deletions < 50 bp (InDels). These variant classes only represent a certain fraction of the genetic variation of an organism. Hence, we reanalyzed whole-genome sequencing data from two experimental populations, which have been divergently selected for FP behavior for over more than 15 generations, performed variant calling for structural variants (SVs) as well as tandem repeats (TRs), and jointly analyzed the data with SNPs and InDels. Genotype imputation and subsequent genome-wide association studies, in combination with expression quantitative trait loci analysis, led to the discovery of multiple variants influencing the GABAergic system. These include a significantly associated TR downstream of the GABA receptor subunit beta-3 (GABRB3) gene, two microRNAs targeting several GABA receptor genes, and dystrophin (DMD), a direct regulator of GABA receptor clustering. Furthermore, we found the transcription factor ETV1 to be associated with the differential expression of 23 genes, which points toward a role of ETV1, together with SMAD4 and KLF14, in the disturbed neurodevelopment of high-feather pecking chickens.

Genetic analysis of production traits and body size measurements and their relationships with metabolic diseases in German Holstein cattle.

This study sheds light on the genetic complexity and interplay of production, body size, and metabolic health in dairy cattle. Phenotypes for body size-related traits from conformation classification (130,166 animals) and production (101,562 animals) of primiparous German Holstein cows were available. Additionally, 21,992, 16,641, and 7,096 animals were from herds with recordings of the metabolic diseases ketosis, displaced abomasum, and milk fever in first, second, and third lactation. Moreover, all animals were genotyped. Heritabilities of traits and genetic correlations between all traits were estimated and GWAS were performed. Heritability was between 0.240 and 0.333 for production and between 0.149 and 0.368 for body size traits. Metabolic diseases were lowly heritable, with estimates ranging from 0.011 to 0.029 in primiparous cows, from 0.008 to 0.031 in second lactation, and from 0.037 to 0.052 in third lactation. Production was found to have negative genetic correlations with body condition score (BCS; -0.279 to -0.343) and udder depth (-0.348 to -0.419). Positive correlations were observed for production and body depth (0.138-0.228), dairy character (DCH) (0.334-0.422), and stature (STAT) (0.084-0.158). In first parity cows, metabolic disease traits were unfavorably correlated with production, with genetic correlations varying from 0.111 to 0.224, implying that higher yielding cows have more metabolic problems. Genetic correlations of disease traits in second and third lactation with production in primiparous cows were low to moderate and in most cases unfavorable. While BCS was negatively correlated with metabolic diseases (-0.255 to -0.470), positive correlations were found between disease traits and DCH (0.269-0.469) as well as STAT (0.172-0.242). Thus, the results indicate that larger and sharper animals with low BCS are more susceptible to metabolic disorders. Genome-wide association studies revealed several significantly associated SNPs for production and conformation traits, confirming previous findings from literature. Moreover, for production and conformation traits, shared significant signals on Bos taurus autosome (BTA) 5 (88.36 Mb) and BTA 6 (86.40 to 87.27 Mb) were found, implying pleiotropy. Additionally, significant SNPs were observed for metabolic diseases on BTA 3, 10, 14, 17, and 26 in first lactation and on BTA 2, 6, 8, 17, and 23 in third lactation. Overall, this study provides important insights into the genetic basis and interrelations of relevant traits in today's Holstein cattle breeding programs, and findings may help to improve selection decisions.

Resistance to chicken amyloid arthropathy is associated with a dysfunctional mutation in serum amyloid A.

Chicken amyloid arthropathy is a debilitating disease with a major impact on animal welfare. Since the disease is triggered by bacterial infection, preventative treatment also contributes to the widespread overuse of antibiotics. Bacterial infection initiates an acute phase response including increased serum amyloid A (SAA) production by the liver. SAA accumulates at sites of infection and in particular in large joints of affected birds. Interestingly, white egg-laying chickens (WL) are resistant to the disease whilst brown egg-laying chickens (BL) are most affected. Disease susceptibility has an immunological basis but the possible contribution of underlying genetic risk factors is not understood. Using a whole genome sequencing approach, we discovered a novel variant in the SAA gene in WL, which is predicted to result in an arginine to serine substitution at position 90 (SAA.R90S). Surprisingly, when overexpressed in chicken hepatocellular carcinoma cells, SAA.R90S was expressed at a higher rate and secreted to a greater degree than the wild-type SAA protein. Moreover, RNASeq analysis showed that the R90S mutant exerted a differential effect on the expression of core transcription factors linked to cell fate determination and cell differentiation. Comparative analysis of gene expression in murine CD4 T-cells stimulated with IL-6/SAA, suggests that SAA.R90S might block an induced cell fate change toward pro-inflammatory T helper 17 cells, which are required for immunological protection against pathogenic bacteria during an acute phase response. Our results provide first mechanistic insights into the genetic resistance of WL to amyloid arthropathy and could be applied to commercial layer breeding programs to improve animal welfare and reduce the negative effects of the overuse of antibiotics.

A multiplexed parallel reaction monitoring assay to monitor bovine pregnancy-associated glycoproteins throughout pregnancy and after gestation.

Bovine pregnancy-associated glycoproteins (boPAGs) are extensively glycosylated secretory proteins of trophoblast cells. Roughly 20 different boPAG members are known but their distribution patterns and degree of glycosylation during pregnancy are not well characterized. The objective of the present study was the development of a parallel reaction monitoring-based assay for the profiling of different boPAGs during pregnancy and after gestation. Furthermore, we investigated the effects of N-glycosylation on our analytical results. BoPAGs were purified from cotyledons of four different pregnancy stages. The assay detects 25 proteotypic peptides from 18 boPAGs in a single run. The highest abundances were found for boPAG 1 in both, glycosylated and deglycosylated samples. Strongest effects of glycosylation were detected during mid and late pregnancy as well as in afterbirth samples. Furthermore, we identified different boPAG-clusters based on the observed relative protein abundances between glycosylated and deglycosylated samples. A linkage between the impact of glycosylation and potential N-glycosylation sites or phylogenetic relation was not detected. In conclusion, the newly developed parallel reaction monitoring-based assay enables for the first time a comprehensive semi-quantitative profiling of 18 different boPAGs during pregnancy and post-partum on protein level, thereby investigating the influence of glycosylation. The results of this study provide new and important starting points to address further research on boPAGs to better understand their physiological role during pregnancy and for the development of new pregnancy detection tests.

Characteristics of thermal images of the mammary gland and of performance in sows differing in health status and parity.

Precision livestock farming can combine sensors and complex data to provide a simple score of meaningful productivity, pig welfare, and farm sustainability, which are the main drivers of modern pig production. Examples include using infrared thermography to monitor the temperature of sows to detect the early stages of the disease. To take account of these drivers, we assigned 697 hybrid (BHZP db. Viktoria) sows to four parity groups. In addition, by pooling clinical findings from every sow and their piglets, sows were classified into three groups for the annotation: healthy, clinically suspicious, and diseased. Besides, the udder was thermographed, and performance data were documented. Results showed that the piglets of diseased sows with eighth or higher parity had the lowest daily weight gain [healthy; 192 g ± 31.2, clinically suspicious; 191 g ± 31.3, diseased; 148 g ± 50.3 (p < 0.05)] and the highest number of stillborn piglets (healthy; 2.2 ± 2.39, clinically suspicious; 2.0 ± 1.62, diseased; 3.91 ± 4.93). Moreover, all diseased sows showed higher maximal skin temperatures by infrared thermography of the udder (p < 0.05). Thus, thermography coupled with Artificial Intelligence (AI) systems can help identify and orient the diagnosis of symptomatic animals to prompt adequate reaction at the earliest time.

eQTL analysis of laying hens divergently selected for feather pecking identifies KLF14 as a potential key regulator for this behavioral disorder.

Feather pecking in chickens is a damaging behavior, seriously impacting animal welfare and leading to economic losses. Feather pecking is a complex trait, which is partly under genetic control. Different hypotheses have been proposed to explain the etiology of feather pecking and notably, several studies have identified similarities between feather pecking and human mental disorders such as obsessive-compulsive disorder and schizophrenia. This study uses transcriptomic and phenotypic data from 167 chickens to map expression quantitative trait loci and to identify regulatory genes with a significant effect on this behavioral disorder using an association weight matrix approach. From 70 of the analyzed differentially expressed genes, 11,790 genome wide significantly associated variants were detected, of which 23 showed multiple associations (≥15). These were located in proximity to a number of genes, which are transcription regulators involved in chromatin binding, nucleic acid metabolism, protein translation and putative regulatory RNAs. The association weight matrix identified 36 genes and the two transcription factors: SP6 (synonym: KLF14) and ENSGALG00000042129 (synonym: CHTOP) as the most significant, with an enrichment of KLF14 binding sites being detectable in 40 differentially expressed genes. This indicates that differential expression between animals showing high and low levels of feather pecking was significantly associated with a genetic variant in proximity to KLF14. This multiallelic variant was located 652 bp downstream of KLF14 and is a deletion of 1-3 bp. We propose that a deletion downstream of the transcription factor KLF14 has a negative impact on the level of T cells in the developing brain of high feather pecking chickens, which leads to developmental and behavioral abnormalities. The lack of CD4 T cells and gamma-Aminobutyric acid (GABA) receptors are important factors for the increased propensity of laying hens to perform feather pecking. As such, KLF14 is a clear candidate regulator for the expression of genes involved in the pathogenic development. By further elucidating the regulatory pathways involved in feather pecking we hope to take significant steps forward in explaining and understanding other mental disorders, not just in chickens.

Structural variants and tandem repeats in the founder individuals of four F2 pig crosses and implications to F2 GWAS results.

BACKGROUND: Structural variants and tandem repeats are relevant sources of genomic variation that are not routinely analyzed in genome wide association studies mainly due to challenging identification and genotyping. Here, we profiled these variants via state-of-the-art strategies in the founder animals of four F2 pig crosses using whole-genome sequence data (20x coverage). The variants were compared at a founder level with the commonly screened SNPs and small indels. At the F2 level, we carried out an association study using imputed structural variants and tandem repeats with four growth and carcass traits followed by a comparison with a previously conducted SNPs and small indels based association study. RESULTS: A total of 13,201 high confidence structural variants and 103,730 polymorphic tandem repeats (with a repeat length of 2-20 bp) were profiled in the founders. We observed a moderate to high (r from 0.48 to 0.57) level of co-localization between SNPs or small indels and structural variants or tandem repeats. In the association step 56.56% of the significant variants were not in high LD with significantly associated SNPs and small indels identified for the same traits in the earlier study and thus presumably not tagged in case of a standard association study. For the four growth and carcass traits investigated, many of the already proposed candidate genes in our previous studies were confirmed and additional ones were identified. Interestingly, a common pattern on how structural variants or tandem repeats regulate the phenotypic traits emerged. Many of the significant variants were embedded or nearby long non-coding RNAs drawing attention to their functional importance. Through which specific mechanisms the identified long non-coding RNAs and their associated structural variants or tandem repeats contribute to quantitative trait variation will need further investigation. CONCLUSIONS: The current study provides insights into the characteristics of structural variants and tandem repeats and their role in association studies. A systematic incorporation of these variants into genome wide association studies is advised. While not of immediate interest for genomic prediction purposes, this will be particularly beneficial for elucidating biological mechanisms driving the complex trait variation.

Replacement of microsatellite markers by imputed medium-density SNP arrays for parentage control in German warmblood horses.

In horses, parentage control is currently performed based on an internationally standardized panel of 17 microsatellite (MS) markers comprising 12 mandatory and five optional markers. Unlike MS, single nucleotide polymorphism (SNP) profiles support a wider portfolio of genomic applications, including parentage control. A transition to SNP-based parentage control is favorable, but requires additional efforts for ensuring generation-overlapping availability of marker genotypes of the same type. To avoid double genotyping of either parents or offspring for changing to SNP technology and enable efficient transition, we tested whether MS genotypes used for parentage control could be reliably imputed from a medium-density SNP panel in German warmblood horses. Imputation accuracy was tested in a tenfold cross-validation with two approaches: within breed (option A) and across breeds (option B). Average imputation accuracies of 97.98% (A) and 96.17% (B) were achieved, respectively. Due to interbreed differences in genotyping rates, five MS markers of low genotyping rate (GTR; < 90%) could be imputed with higher accuracy within breed (98.18%) than across breeds (90.73%). MS markers with high GTR performed homogeneously well in option B (98.44%) and showed slightly lower accuracy in option A (97.90%). Among these markers, AHT5 proved to be problematic for imputation regardless of the approach, revealing accuracies of 86.40% (A) and 88.70% (B). Better results for MS markers with high GTR and savings in computational processing justified the choice of option B for routine implementation. To date, more than 9500 horses have undergone the new parentage control based on imputed MS genotypes.